Fast growing double tearing modes in a tokamak plasma 
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Configurations with nearby multiple resonant surfaces have broad spectra of linearly unstable 
coupled tearing modes with dominant high poloidal mode numbers m. This was recently shown for 
the case of multiple q = 1 resonances [Bierwage et al., Phys. Rev. Lett. 94 (6), 6500f (2005)]. In 
the present work, similar behavior is found for double tearing modes (DTM) on resonant surfaces 
with q > 1. A detailed analysis of linear instability characteristics of DTMs with various mode 
numbers m is performed using numerical simulations. The mode structures and dispersion relations 
for linearly unstable modes are calculated. Comparisons between low- and higher-m modes are 
carried out, and the roles of the inter-resonance distance and of the magnetic Reynolds number 
5h p are investigated. High-m modes are found to be destabilized when the distance between the 
resonant surfaces is small. They dominate over low-m modes in a wide range of Sh p , including 
regimes relevant for tokamak operation. These results may be readily applied to configurations with 
more than two resonant surfaces. 
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I. INTRODUCTION 

Double tearing modes (DTM) arc coupled tearing 
modes on adjacent resonant surfaces which effectively 
"drive each other" Q, 0, Q. The resulting instability 
is stronger than a single tearing mode (STM), whereby 
modes with high poloidal mode numbers m may become 
dominant Q. 

DTMs are thought to be involved in a variety of dy- 
namical processes in tokamak plasmas. In the current- 
ramp-up phase of the tokamak operation they were used 
to explain the anomalously strong current penetration 
ULEIlI an d short-wavelength MHD activity (e.g., Mirnov 
oscillations) @, Experimental observations of par- 
tial and compound sawtooth crashes (internal disrup- 
tions) n us m m eh may be explained through 

kink-tea ring modes on double/multiple resonant surfaces 
[Til [TEL Il6l Il7j. Phenomena associated with low-beta 
disruptions in the presence of double resonant surfaces 
with q s > 1 may also be understood in terms of DTM 
activity [3, EH O In some way related to DTMs is 
the coupling between modes on resonant surfaces with 
different q s , which may explain some salient features of 
major disruptions 0, l20l . l2lj . A more recent application 
is related to the formation of internal transport barriers 
(ITB) in advanced tokamaks. DTMs were suggested to 
be involved in this process due to the fact that ITBs arc 
often observed in the vicinity of resonant surfaces or near 

Qram 



the distance between the resonant surfaces is sufficiently 
small @- This is reflected by the linear growth rate 
7iin scaling with the magnetic Reynolds number Sh_ p as 



It was found that DTMs with any mode number m 
become similar to m — I internal kink modes l25l when 
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7ii n oc ■ When the distance between the resonant 

surfaces is large, the structure of the instability resem- 
bles that of individual STMs localized on each resonant 
surface. In this case, the scaling law is jn n oc 5" H p^ 5 
[26| . Furthermore, it was found that the coupling in a 
DTM may also destabilize linearly stable tearing modes 
EH and speed up the growth of magnetic islands in 
the nonlinear regime El HI The instability of 
DTMs may be enhanced by effects like anomalous elec- 
tron viscosity 0, |2^, finite beta and sheared toroidal 
flows j^ij. Poloidal shear flows, on the other hand, may 
have a stabilizing effect [Tel l3fl l3^ |. raising the possi- 
bility of a dynamical interplay between shear flows and 
DTMs [Til l33l | . This idea is supported by recent studies 
which indicate that MHD activity associated with multi- 
ple resonant surfaces may produce a significant amount 
of sheared poloidal flows [23 . l3(| . 

Until now, studies related to DTMs in cylindrical (or 
toroidal) plasmas were focussing mainly on the role of 
the modes with the lowest poloidal and toroidal mode 
numbers, m and n. High-m tearing modes were associ- 
ated only with correspondingly high values of q s = m/n 
HE HQ. Recently it was demonstrated for the case of 
q s = 1 double and triple tearing modes that modes with 
high m may dominate over low-m modes in a cylindrical 
plasma even on low-q resonant surfaces Q . There it was 
found that in configurations where the distance between 
neighboring resonant surfaces is sufficiently small, broad 
spectra of unstable modes exist and that the dominant 
modes can have m ~ 0(10). 

Motivated by the results of Ref. Q , in the present work 
we investigate the conditions under which broad spectra 
of multiple tearing modes with dominant highcr-m modes 
may arise. For this purpose numerical simulations based 
on the reduced resistive magnctohydrodynamic model in 
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cylindrical geometry are used. The dependence of the lin- 
ear growth rate on parameters such as the inter-resonance 
distance, the dissipation coefficients as well as the mag- 
netic shear will be addressed. 

The main results of the present work are the following. 
It is shown that broad spectra of unstable DTMs are also 
found in configurations with q s > 1, thus generalizing the 
results of Ref . . The instability of modes with m > 1 
depends strongly on the distance between adjacent res- 
onant surfaces. For the resistivity dependence, a scaling 
law 7ii n oc is found with 1/3 < a < 3/5, which is 
valid in a range of high values of Sn p . This is in agree- 
ment with former studies 0, l28j . For intermediate values 
of S^p n o power law is identified and for low Sjj p it is 
found that 7i; n oc Sn p independently of m. For high m 
and small D\2 the range where the relation 7i; n oc 
does not hold can be rather wide, and it is especially in 
this regime where high-m modes are found to be dom- 
inant. This regime may include ranges of Sh p that are 
relevant to tokamak operation. 

The present paper is organized as follows. Section [H] 
contains an introduction of the model equations and the 
numerical method used. Numerical results are presented 
in Section II I II followed by a discussion and conclusions 
in Section llVl 



II. MODEL 
A. Reduced MHD equations 

The reduced set of magnetohydrodynamic (RMHD) 
equations in cylindrical geometry [35| is used. It is ob- 
tained from the full MHD model through a high-aspect 
ratio and low-/? expansion: e = a/Ro <C 1, <C 1 (e.g., 
Here, a and Ro are, respectively, the minor and 
major radius of the torus, e is the inverse aspect ratio, 
and the plasma beta (3 — p/(Bq/2[1o), is the ratio of 
thermodynamic pressure to magnetic pressure. 

The high-aspect-ratio ordering and the presence of a 
strong axial magnetic field Bq allow to express the mag- 
netic field in terms of a magnetic flux function ^ as 
B = Bqz + x z. Here, z is the unit vector in the axial 
direction and Bo is taken to be constant. In the MHD 
ordering, the single-fluid velocity V is approximated by 
the E x B drift velocity V E = — V$ x z/B a , where $ 
is the electrostatic potential (or stream function) . In the 
limit of zero (3 the pressure equation is decoupled and the 
RMHD equations take the form of a two-field model: 

= -B V$ + — Vi^ + E z , (1) 
Mo 

p m d t vi$ = -— b- vvi* + i/ m vivi$. (2) 

Mo 

Here p m is the mass density (constant due to the assump- 
tion of incomprcssibility) , r\ the plasma resistivity and v m 
is the kinematic viscosity. The time-independent electric 



field E z satisfies VE Z x z = 0. The convective derivative 
is defined as d t = d t + V E • V. 

Cylindrical geometry is chosen, which gives the right- 
handed set of coordinates (r, z) where r G [0, a] is the 
radius, $ G [0, 27r] the poloidal angle and z £ [0, 2nR ] 
the axial coordinate (related to the toroidal angle if via 
z = —Rocp). The Poisson bracket is defined as [f,g] = 
-(d r fd$g — d r gd$f), and the Laplacian is approximated 
by V 2 ~ Vj_ = \d r rd T + In this ordering, toroidal 

topology is preserved by retaining periodicity in the axial 
coordinate z and the poloidal angle &. 

Normalizing the time by the poloidal Alfven time 
t~h p = y/fiap m a/ B$(a), the radial coordinate by a, 
and introducing the angular coordinate £ = z/R q(a) 1 
the normalized variables are ip = \& /aB^(a) and <f> = 
<I>/(a 2 /THp). Note that f is related to the toroidal angle 
ip via C = —tp/q(a), with q(a) = aB^/ RqB§{cl) being the 
safety factor at r = a [c.f., Eq. JSJ]. The normalization 
for the source term is Eq = E z /\rjoB^{a) / fioa}. With 
these normalizations the Eqs. JTJ and J5J become 



d t u 



'Hp 

d c3 + Re^lV 2 , 



(3) 
(4) 



where the current density j and the vorticity u are re- 



lated to i(j and (j) through j 



-Vi ip and u = 



respectively. The magnetic Reynolds number S'jjp (also 
called Lundquist number) and the kinematic Reynolds 
number i?ejj P arc defined as S^p = t v /th p and Ren p = 
a 2 / (v m TR P / Pm) , respectively. Here, t v = a 2 po/r]{ r = 0) 
is the resistive diffusion time. The resistivity profile 
is given by fj(r), normalized such that fj(r = 0) = 1. 
The relative strengths of kinetic viscosity and (resis- 
tive) diffusion are characterized by the Prandtl number 
Pr = Sjjp/^ejjp oc v/r/, with v = v m / p m being the spe- 
cific ion viscosity. 



B. RMHD equilibrium and resonant modes 

The equilibrium state is defined as dtip = dt4> = 0. 
Since the effect of the plasma pressure is neglected, the 
structure of the equilibrium magnetic field B cq is given 
in terms of the tokamak safety factor q, defined as 



q 



B eq ■ VC 

B oq - W 



(5) 



The safety factor measures the field line pitch by counting 
how many times a magnetic field line goes the long way 
around the torus (27ri?o) after one turn the short way 
around (2wr). Here, q = q(r), so that the equilibrium 
magnetic flux surfaces, where ip = const., are uniquely 
defined by the radius r. Each field variable / is written 
in terms of a time-independent equilibrium component / 
and a time-dependent perturbation / as 



(6) 
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It is assumed that the equilibrium state is free of flows. 

= u = O. (7) 

In general, a tokamak plasma has magnetic surfaces 
where q s = q(r s ) = m/n, with integers m and n. These 
are called rational or resonant surfaces. The radius r s is 
called resonant radius, because an infinitcsimally small 
helical magnetic perturbation with hclicity h — m/n is 
resonant with the magnetic field structure in the vicin- 
ity of r s . Such a resonant perturbation 5ip miTl does not 
bend field lines, so there is no restoring force. If the per- 
turbation leads the system to a state of lower energy, the 
amplitude of the resonant perturbation will grow and the 
mode is said to be unstable. 



C. Fourier representation 

In order to study the properties of such resonant modes 
it is useful to apply a Fourier transform with respect to 
the periodic coordinates: (•&, — > (m, n). This repre- 
sentation also gives an efficient and accurate numerical 
model. Substituting for each field variable / in Eqs. © 
and the Fourier expansion 



f(r, 0, C, t) = \ fm,n( r > t)-e l(m "-" c) + 



ex., 



(8) 



one obtains equations for the individual Fourier modes, 



'-'Hp (vjrii,n Fm,n) i 

d t u m , n = [u, 4>] m n + [j, tp] m n - mj m ,„ 



(9) 
(10) 



where V 2 



^d r rd r — m 2 /r 2 and the nonlinear terms 
[f,g] m ,n have acquired the form of convolutions. In this 
model, if a perturbation is applied only to modes of a 
given hclicity h = m/n, modes with different helicitics 
will not be excited and the problem is effectively reduced 
to a two-dimensional one. In the following, we will ex- 
clusively refer to individual Fourier components of the 
field variables and usually omit the (m, n) subscripts for 
convenience. Note that the equilibrium fields / = f(r) 
have only (m, n) = (0, 0) components. 



D. Equilibrium model 

In order to study configurations with two resonant sur- 
faces, the following model formula for the equilibrium q 
profile is used: 



l(r)=q .F 1 (r).{l + (r/r ) 2 ^} 



(11) 



FIG. 1: (Color online) Safety factor profiles with two resonant 
surfaces used in this study. The Case (la), plotted in (a), 
has two q s = 1 resonances, located at r = r s i and r S 2, a 
relatively small distance D12 = 0.06 apart. In (b), variants 
of this profile, with different inter-resonance distances D12 
and different q B are shown. The parameter values required to 
reproduce Case (la) are given in Table Q and the geometric 
characteristics of all cases are listed in Table ITT1 



where 



'0 
u(r) 
Fi(r) 



[m/{nq )] 
.2 



A»(»"a) 



-l/2/i(r A ) 



1 + h exp j - [(r - rn)/ri 2 ] 2 | 



The parameter set {<7o, r Ai MOi n i m i n } is used to de- 
sign the underlying monotonic profile. The parame- 
ters {/i,T , n,ri2} describe the Gaussian bump which is 
used to create non-monotonic profiles with two (or three) 
q s = m/n resonant surfaces. The DTM configurations 
used in this study are shown in Fig.^ The corresponding 
model parameters are listed in Tableland the geometric 
characteristics of all cases are given in Table |n] 

With q given by Eq. (|ll|l . ip and j were calculated 
using the relations 



" X = --? and 
r dr 



(12) 



The resonant surfaces where q = q s are labeled with 
r s i and r S 2, and the inter-resonance distance is given by 
D12 = \r S 2 — r s i\. The magnetic shear profile is defined 

as 



r dq 
q dr 



dflng) 
d(lnr)' 



(13) 



and the local magnetic shear at the resonant radius r S i is 
&i s(rg^). 

The time- independent source term compensates the 
resistive dissipation of the equilibrium current profile, 
i.e., Eq = rjj (diffusive equilibrium). For numerical sim- 
ulations where the temporal evolution of the resistivity 
profile r] is neglected one often assumes fj(r) = j(r = 
0)/j(r). This will generally lead to different values of 
the resistivity at different r si . In order to simplify the 
comparison between growth rates of modes associated 
with different resonant surfaces, a homogeneous resistiv- 
ity profile, f)(r) — 1, is used in this study. Numerical 
tests indicate that the details of i) have no significant ef- 
fect on any of the qualitative characteristics discussed in 
this paper. 



E. Linearized equations 

When the amplitudes and gradients of the perturbed 
fields are sufficiently small the nonlinear terms in Eqs. @ 
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go r A Mo in m n /i m n 2 

1.3 0.655 3.8824 1 1 -0.238 0.4286 0.304 



TABLE I: Parameter values for Eq. 1111 . giving the q profile 
shown in Fig. ^a). The other cases are readily obtained by 
changing m and n (so that q B = m/n) and go- 



Case 


g s 


Qmin 


D 12 


Sl 


S2 


(la) 


l 


0.99 


0.06 


-0.10 


0.12 


(lb) 


l 


0.96 


0.21 


-0.20 


0.45 


(Ic) 


l 


0.92 


0.31 


-0.20 


0.66 


(II) 


3/2 


1.49 


0.06 


-0.10 


0.12 


(Ilia) 


2 


1.92 


0.31 


-0.20 


0.66 


(nib) 


2 


1.99 


0.06 


-0.10 


0.12 


(IV) 


5/2 


2.49 


0.06 


-0.10 


0.12 


(V) 


3 


2.99 


0.06 


-0.10 


0.12 



FIG. 2: (Color online) Numerical convergence of the linear 
growth rates 7ii n (m) with increasing N r (number of radial grid 
points per unit length). The configuration used is Case (la) 
with D 12 = 0.06, S Hp = 10 s and i?e Hp = 10 12 . The growth 
rates are shown for modes m — 1, 2, 5 and 8. The m = 5 mode 
is dominant in this case. The small symbols represent EVP 
results, calculated in the range 250 < N r < 1430. The large 
symbols represent IVP result. The IVP solver was run with 
up to N r = 3000. The agreement between both methods and 
the numerical convergence is clearly shown. The horizontal 
dotted lines indicate the linear growth rate values for N r = 
730. This is the standard iV r -value used for Sh p = 10 8 in this 
paper (cf., Table ITTTl . It can be seen that these growth rates 
differ from the N r = 3000-results only by a few percent. This 
accuracy is sufficient for the purpose of this paper. 



5h p : 


< io 5 


10 5 -» 10 7 


10 7 -» 10 8 


> 10 8 


N r : 


180-360 


430-570 


730 


1460 



TABLE II: Geometric properties of the q profiles shown in 

Fig. Eft) . 



and HI Up may be neglected and one obtains the linearized 
RMHD equations. In the linear system, the time depen- 
dence of a perturbed field variable / takes the form 



/(r,t) =/(r).exp(Ai), 



(14) 



where A is a complex number. Using Eq. I|14|) . the system 
of equations © and l|10[l becomes 



A</> 



Xu = i yri — —j ]^-d r rd r — —j-j ip 
im ( sis — 2) d r s\ ~ 

+—[- L - — )i> 

r \ rq q ) 

i n a a m 2 

— OrTOr 7T I U. 



(16) 



Re-Hp 

Equations l|15|) and (|16|l are obtained by applying 
Eqs. (|12f) and and expressing the equilibrium fields 
in terms of the safety factor q = q(r) and the magnetic 
shear s = s(r). 

Modes for which the linear growth rate is positive, 
7iin = 5i{A} > 0, are said to be linearly unstable. Their 
amplitudes grow exponentially in time, e.g., 



(17) 



(in RMHD, 3{A} = 0). The linear growth rate 7n n is a 
function of the mode numbers. 



7lin = 7lin(™,«) 



(18) 



(spectrum of growth rates), and written this way it is 
generally referred to as the dispersion relation. Since this 



TABLE III: Typical values for N r (number of radial grid 
points per unit length) used in the EVP solver. The IVP 
code is usually run with 2000 < N r < 3000. 



study is restricted to modes with unique hclicity, h, it is 
sufficient to specify m for a given q s = h (so n = m/q s ). 
Along with the radial structure of the cigenmodcs, ip{r) 
and 4>{r), the linear growth rates are our main tool for 
characterizing the linear instability of DTMs under vari- 
ous conditions. 



F. Numerical method 

After discretizing Eqs. I|15|) and l|16|) with respect to 
the radial coordinate, the linearized RMHD model may 
be written as a generalized eigenvalue problem, 



BAx = Ax, 



(19) 



where x 



A is a diagonal matrix containing 

the eigenvalues A, and A and B are coefficient matri- 
ces, which also include the finite-difference operators. All 
eigenmodes Xj with corresponding eigenvalues Xj are ob- 
tained using an eigenvalue problem (EVP) solver. The 
accuracy of the results is checked by running tests with 
different numbers of radial grid points and through com- 
parison with results obtained by solving the linearized 
version of Eqs. JTJ and (J2J as an initial value problem 
(IVP). The numerical convergence and the agreement 
between the EVP and IVP results are demonstrated in 
Fig- El I 11 Tablc lTTll the radial resolution used for different 
values of Su p is specified. Due to numerical constraints, 
the regime 5h p > 10 10 was not accessible with sufficient 
accuracy. 
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FIG. 3: (Color online) Mode structures of the unstable q s = 
1 DTM eigenmodes of Case (la), where D 12 = 0.06. The 
locations of the resonant surfaces are indicated by vertical 
dashed lines. Diagrams (a) and (c) show, respectively, if) and 
4>/r of the m = 1 mode. In (b) and (d) the mode structures 
for m — 8 are shown, which are representative for other modes 
with m > 1. In our notation, the eigenmode M^ 1 ' extends to 
r = r s i, and M' 2 ' to r S 2. These results are obtained with 
5h p = 10 6 , but similar mode structures are also found for 
other values of Sh p > 10 5 . The linear growth rates 7^ and 
7un are a l so shown in (b) and (d). 



FIG. 4: (Color online) Mode structures of the unstable q s = 
2 DTM eigenmodes with (m,n) = (2,1). In (a) and (c), 
respectively, 1]} and (j>/r are shown for Case (Ilia) (D12 = 
0.31). In (b) and (d) the results for Case (Illb) (D 12 = 0.06) 
are plotted. These results were obtained with Sh p = 10 6 . 

III. RESULTS 

A. Mode structures 

The radial structure of an eigenmode determines which 
part of the plasma is affected by the instability, since 4>(r) 
is related to the radial and poloidal components of the 
plasma displacement velocity, 

v r = cf> and v.ff = d r (f>. (20) 

r 

In general, for each (to, n) there may be as many unsta- 
ble eigenmodes as there are resonant surfaces q s = m,/n. 
The radial mode structures of unstable q s = 1 DTMs are 
shown in Fig. |3J The configuration used is Case (la), 
where D 12 = 0.06. For m = 1 [Fig. EJa) and (c)] there 
are two unstable eigenmodes. For m > 1 [Fig. |3fb) and 
(d)] only one eigenmode is unstable (the to = 8 mode 
structure shown is representative for other modes with 
m > 1). In our notation, the eigenmode is associ- 

ated with the innermost resonant surface r s \. For m = 1 
it has a finite amplitude in the region < r < r sl . The 
eigenmode is active in the region < r < r S 2 for 

m = 1, and mainly in the region r sl < < r S 2 for to > 1. 

In Figure 01 the mode structures of q s = 2 DTMs 
with (m, n) = (2, 1) are plotted. Results are shown for 
two cases with different distances Dyi'. Case (Ilia) with 
D 12 = 0.31 [Fig. Ufa) and (c)], and Case (Illb) with 
D12 = 0.06 [Fig. 0Jb) and (d)]. In agreement with pre- 
vious works on DTMs, it is found that two individual 
modes are present in the case of larger D12: M^ 1 ' with 
even parity around r sl and with odd parity psl l37l| . 
In the case of smaller D12 only the eigenmode M&> is 
found to be unstable. 

The eigenmode is essentially an STM associated 
with the r = r sl resonant surface (here, a negative-shear 
surface) , since it is practically unaffected by the presence 
of resonant surfaces beyond r = r s i. The actual DTMs 



FIG. 5: (Color online) Dependence of the linear DTM growth 
rates on the Prandtl number [Case (la)]. The growth rates of 
eigenmodes with m — 1, 2, 8 are plotted for fixed Sh p = 10 6 . 



FIG. 6: (Color online) Comparison between dispersion rela- 
tions for (a) q s — 1, 2 and 3, and (b) g s = 3/2 and 5/2. For 
all cases D12 = 0.06 (cf., Table ITTt . For m — 1, circles in 
diagram (a), growth rates of both unstable modes are shown: 

Tito = 4A x 10 ~ 3 ' Tiin = 16 x 10 ~ 3 - A11 other growth rates 
belong to M^'-type eigenmodes, since the M^'-type modes 
are stable for m > 1 in these cases. For all cases Sh p = 10 6 . 

(for sufficiently small D12) are the M^> eigenmodes. The 
radial structure of DTMs with to > 1 is very similar for 
different q s , as may be seen by comparing the profiles (b) 
and (d) in Figs. and 0] Let us note that Cases (II), 
(IV) and (V) of Table El where q s = 3/2, 5/2 and 3, 
respectively, have eigenmode structures very similar to 
those found for q s = 2 in Fig. 

The local resistivity rj(r S i) and the magnetic shear 
s(r s i) at a resonant surface r S i, together with the dis- 
tance between neighboring resonant surfaces -D12, deter- 
mine which eigenmode, M' 1 ' or M^ 2 \ will be the domi- 
nant mode for a given m. An exception is the STM-likc 
eigenmode M^- 1 ': as its mode structure indicates, it is 
not affected by D12. 



B. Viscosity effect 

Although, it is difficult to obtain reliable values for the 
ion viscosity in the tokamak core, it is usually thought 
to be very low. However, it is required in most nonlin- 
ear simulations for the purpose of providing a cut-off at 
short wavelengths in order to be consistent with the fi- 
nite number of grid points or Fourier modes. As may be 
seen from Fig. [S] as long as the Prandtl number satisfies 
Pr < 0.1, the viscosity has practically no effect on the 
linear growth rates of DTMs. A similar result is also ob- 
tained for other values of Sh p as well as for single and 
triple tearing modes, so it may be considered a generic 
characteristic of tearing modes. The mode structures in 
Figs. |3 and 0| and all the following results were obtained 
in the regime where Pr < 0.01, so the viscosity effect 
will not be discussed further in this paper. Let us note 
that the stabilizing effect of viscosity for Pr > 1, evident 
in Fig. [31 was pointed out previously by Ofman |3lj | in 
a study on DTMs in the presence of equilibrium shear 
flows. 



C. Dispersion relations for q s > 1 

Recently it was discovered that configurations with 
nearby q s — 1 resonant surfaces have broad spectra of 
linearly unstable modes Q. In particular, if the distance 
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FIG. 7: (Color online) Role of D12 (variable shears) for the 
instability of q s — 1 DTMs. In (a) the q profiles are shown 
[Cases (Ia,b,c) in Table ITT1 and in (b) the corresponding dis- 
persion relations are plotted. For all cases Sh p = 10 6 , and 
only the growth rate of the dominant mode is shown for each 



FIG. 8: (Color online) Role of D\2 (constant shears) for the 
instability of q e = 1 DTMs. In (a) the q profiles used are 
shown. The basic profile is Case (lb), with D12 = 0.21. The 
other profiles (D12 < 0.17) were obtained by gradually reduc- 
ing D12, but holding the shears Si = —0.20 and S2 = 0.45 
constant. For this, the model equation (II IP was multiplied 
by another factor, -Fb(r), defined equivalently to Fi(r). In 
(b) the corresponding dispersion relations are shown. For all 
cases Sh p = 10 6 , and only the growth rate of the dominant 
mode is shown for each m. 



between the resonant radii is sufficiently small, modes 
with m ~ O(10) were found to have linear growth rates 
several times higher than the m = 1 mode. 

Indeed, similar behavior is observed when the reso- 
nances are located at higher q values, such as q s = 3/2, 
2, 5/2 and 3, corresponding to Cases (lib), (III), (IV) 
and (V), respectively. The dispersion relations for these 
cases are shown in Fig. [fj] All these q profiles have the 
same values for D12, s\ and S2- They only differ from 
each other in the values of q m in and q s . Let us remark 
that for all cases the width of the spectrum as well as 
the location of the maximum growth rate is practically 
identical: 7i; n > for 1 < m < 18, and the peak 



Tpeak = 7iin(m peak ) = Max{7 Un (m)} 



(21) 



is located at m = 8. Except for the m = 1 mode [q s 
1 . Case (la)], all growth rates in Fig. are associated 
with -type eigenmodes. Due to these results it may 
be conjectured that all DTMs [except for M^(m = 1)] 
behave similarly for any q s , including q s = 1. 



D. Role of Di 2 

Consider the scenario where <? m i n gradually drops be- 
low q s = 1 due to an increase in the current density. 
During this process the distance between the resonant 
surfaces will grow in time, starting from D12 = when 
5min = q s - The series of Cases (la), (lb) and (Ic), shown 



in Fig. (a), represent the equilibria at successive in- 
stants in time during such a process. Similar cases were 
previously investigated in Ref. where the effect of 
D12 on 7iin(m = 3,n = 1) and the corresponding scal- 
ing exponent a (in 7u n oc S^) was characterized. In 
Ref. [3^], the role of g m i n was studied in the context of 
partial and full rcconncction. 

The dispersion relations for the three cases of Fig. 0a) 
are plotted in Fig. 0h) . Clearly, the width of the spec- 
trum is reduced as D12 increases, and the m = 1 mode 
eventually becomes dominant. 

Note that Cases (Ia)-(Ic) all have different magnetic 
shears s± and S2 [cf., Table ITTj . The effect of varying 
only D12 is illustrated in Fig. EI Starting with Case (lb), 
where D12 = 0.21, the q profile is gradually modified in 
such a way that D\2 decreases down to 0.08, while both 
shears, s\ and S2, are held constant. The profiles used 
are plotted in Fig.|SJa) and the corresponding dispersion 
relations in Fig. Efb). These results show that D\2 con- 
trols the broadness of the spectrum. Modes with higher 
m and with higher growth rates than the lowcr-m modes 
appear when D12 is decreased. 

It is known that a tearing mode is stable when the 
local magnetic shear at the resonant surface is zero [26| . 
The results in Fig. indicate that for DTMs with m > 1 
the destabilizing effect of small D12 dominates over the 
effect of low local shears [|si| and S2 become small when 
Aq = q s — q min is reduced (cf., Table HT) l]. 

The dependence of the DTM growth rates on the pa- 
rameter q m i n (meaning, simultaneous variation of -D12, 
si and s 2 ) was studied previously by Ishii et al. j2^| for 
the (m, n) = (3, 1) mode and a set of q profiles similar to 
those in Fig. 0a). They found that the curve 7im(9min) is 
not monotonic. In Fig.^a) it is shown that this is also 
true for the dependence of 7ij n on D12 alone (constant 
shears): in the range shown (0.08 < D12 < 0.21), the 



growth rate of the M^(m = 1) mode, 7^ 



(2) 



1), has 



a maximum around D12 = 0.14. A comparison between 
the results for 5 Hp = 10 6 in Fig. Efa) and 5 Hp = 10 7 
in Fig. EJb) shows that the value 5h p influences the lo- 
cation of the peak. This is most likely related to the 
distance D12 becoming comparable to the linear resistive 
layer width S v , which is estimated after Ref. , 



7im(«i) 



(m/r„ 



'Hp 



1/4 



(22) 



FIG. 9: (Color online) D12 dependence of the growth rates of 
m = 1 and m = 2 DTMs with q s = 1. While varying D12, 
the shears si = —0.20 and S2 = 0.45 were held constant. The 
values for Lundquist number are (a) Sn P = 10 6 and (b) 5h p = 
10 7 . For several data points, indicated by arrows, estimates 
are given for the linear resistive layer width S v defined by 
Eq. 



where B' s = s(r s i)/q(r S i). In the following section, the 
role of Sup will be studied in more detail. 

It is noted that the growth rate 7^ \{m = 1) does not 
depend on -D12, as it is expected from the mode structure 
MW(m = 1) [Fig. I^a) and (c)]. Hence, the increase of 
the m = 1 growth rate in Fig. 0b) is caused by the 
increase in the magnetic shear |si|. 
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FIG. 10: (Color online) DTM dispersion relations for Sh p = 
10 6 , 10 7 and 10 8 obtained with Case (la) where q B — 1. 
For m = 1 there are two unstable eigenmodes, and 
M (2) . The growth rates of both the dominant (M (1) ) and 
the secondary (M^) m = 1 eigenmode are shown for 
5h p = 10 6 (they almost coincide for higher Sh p ). Verti- 
cal dashed lines indicate the locations m pca k of the peaks 
7 P cak = Max{7ii n (m)}. Similar results are obtained for cases 
with q s > 1. 



E. Role of resistivity 



FIG. 11: (Color online) Sh p dependence of the linear growth 
rates of the two m = 1 DTM eigenmodes M (1) and M (2) 
of Case (la) [cf., Fig. |3] (a) and (c)]. For comparison the 
growth-rate scalings of "corresponding" STMs are shown as 
well. This means, that the instability of the DTM eigenmode 
Af is compared with the STM mode on the resonant surface 
(r s i,si) (monotonic profile with negative shear), while M*- 2 ' 
is compared with the STM mode on (r S 2,S2) (positive-shear 
profile). The STM results are plotted as straight lines, fitted 
to the actual data. In the regime 5h p ^ 10 6 , the 7n n oc 
scaling with a — 1/3 is followed closely by both STMs: a w 
0.32. 



In the previous sections the linear instability of DTMs 
was investigated at Sh p = 10 6 , since this value lies in the 
regime where nonlinear simulations of MHD instabilities 
arc often performed: 10 4 < Sh p < 10 7 . While large toka- 
maks typically operate in regimes where Sh p <^ 10 8 (ex- 
cept in the early current-ramp- up phase), it is difficult to 
access this low-collisional regime with nonlinear simula- 
tions. Therefore, the linear study of the Sh p dependence 
of DTM growth rates is an important tool for relating 
nonlinear simulation results to plasma conditions. 

In Fig.^|the dispersion relations of Case (la) is shown 
for S Hp = 10 6 , 10 7 and 10 8 . The growth rates of all 
modes are found to decrease with increasing Sh p , as it is 
expected for resistive instabilities. Let us remark that the 
width of the spectrum is not affected by the variation of 
Sh p : unstable modes are found in the range 1 < m < 18. 
However, it can be seen that the growth rates of modes 
with higher m drop more rapidly than growth rates of 
low-m modes, when Sh p is increased, which is also re- 
flected by the shift of the peak 7 pea k to lower m. This 
is important, because the mode number of the fastest- 
growing mode determines the size of the magnetic islands 
formed in the early nonlinear regime Q. 

In the remaining part of this section, the Sh p depen- 
dence of DTM growth rates is examined in detail. Com- 
parisons between low- and higher-m modes are made and 
the role of the distance D\ 2 will also be emphasized. 
First, cases where <j s = 1 are considered, whereby a dis- 
tinction is made between DTMs with m — 1 and m > 1. 
The results obtained for q s = 1 are readily applied to 
cases with q s > 1, as will be shown at the end of this 
section where q s = 2 DTMs are considered. 

DTMs with m = l 

It is known that DTMs with m — 1 behave similarly to 
rn = 1 STMs, regardless of the location of the resonant 
surfaces and their mutual distance [||. Due to this close 
relationship between the m = 1 STMs and the m = 1 
DTMs it is possible to decompose the non-monotonic q 
profile into two monotonic ones — one with a resonance 
at r s i and negative shear s\, the other with a resonance at 
r s2 and positive shear s 2 — and compare the growth rates 
and S'Hp dependences of the STM eigenmodes with the 



FIG. 12: (Color online) Sh p dependence of the linear growth 
rates of the two m — 1 DTM eigenmodes, (a) for D12 = 0.31 
[Case (Ic)] and (b) for D12 = 0.21 [Case (lb)]. Corresponding 
STM data are shown as well, as in Fig. 1111 



corresponding DTMs. Our aim is to characterize devia- 
tions between STMs and DTMs in certain ranges of Sh p - 
For clarity, both AfW- and M( 2 )-type eigenmodes in a 
double-tearing configuration are referred to as "DTMs." 

In Fig. ^2 this comparison is performed for Case (la), 
where D 12 = 0.06 [Fig. □ (a)] . The results for Cases (lb) 
and (Ic), where D i2 = 0.21 and D i3 = 0.31, respectively, 
are presented in Fig. 

In Case (la), |si| w s 2 (cf., Tablc HT)) . so that the growth 
rates of the two STMs are almost equal (Fig. Ill|) . Both 
have a = 0.32 ss 1/3, in agreement with linear theory It 
can be observed that the DTM growth rates coincide with 
the corresponding STM growth rates only at relatively 

high values of Sh p : 7^ f° r ^Hp > 10 7 and 7^ for Sh p > 
10 s . The agreement between m = 1 DTMs and STMs 
for intermediate S'Hp is significantly improved when the 
distance between the resonant surfaces is increased. This 
can be seen by comparing Fig.^^whh Fig. 1121 

The linear growth rate of a tearing mode plotted as 
a function of the magnetic Reynolds number, 7r m (S l H P ), 
always has a maximum at a certain value Syv v = <S max . 
For instance, for 7^ in Fig.^] one finds S mayi «5x 10 5 . 
In the regime where S'Hp < S max the width of the linear 
resistive layer S v [Eq. {25J] is comparable to D\ 2 . Due to 
the strong resistive diffusion in the regime Sh p < S max , 
the two resonant surfaces are effectively seen as a single 
q s = 1 "surface." This leads to the reduction in the 
DTM growth rate apparent in Fig. ^] For Sh p <C S max 
a linear dependence 7i; n oc Sh p (i.e., a = — 1) is obtained 
independently of m. 

In summary, for m = 1 DTMs it is found that for suf- 
ficiently small D12, the growth rate of the M^ 1 mode 
exhibits no scaling law in a wide range of Sh p - For 
Sh p -C S max a linear dependence is observed and for 
Sh p S> S max the scaling exponent a = 1/3 is obtained 
for both eigenmodes. However, the range of Sh p in which 
deviations from the a = 1/3 scaling are observed may ex- 
tend to high Sh p when D12 is small. 
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FIG. 13: (Color online) Su P dependence of the linear growth 
rates of q s — 1 DTMs with m = 1,2,3 for a small distance 
D12 = 0.06 [Case (la)]. The dotted lines represent the scaling 
law 7n n = Sjjp with a = 1/3 and 3/5. 



FIG. 14: (Color online) Sn p dependence of the linear growth 
rates of q s — 1 DTMs with m = 1, 2 for the distance D 12 = 
0.21 [Case (lb)]. 



DTMs with m > 1 

Here the Sh p dependence of q s = 1 DTMs with m > 1 
is analyzed. First, consider Case (la) where the inter- 
resonance distance is small: D12 = 0.06. As indicated 

in Fig. we have selected the growth rates 7 lin (m = 
(2) 

2) and Tun ( m = °), and plotted them as functions of 
5 Hp in Fig. US It is noted that in Case (la) only A/ (2) - 
type modes are unstable for m > 1. For comparison, 

7un( m = 1) ana - 7ihi( m = 1) fr° m Fig. 1111 are shown 
as well. It can be seen that, similarly to the m = 1 

(2) (2} 

mode, 7j in ; (m = 2) and 7 lin (wi. = 8) have a = — 1 for 
<Sh p < -Smax, where S max (m = 2) w 10 5 and S max (m = 
8) w 2 x 10 5 . In the range of S*h p shown in Fig. IPjI 

(2) 

only 7 lin (m = 8) approaches the scaling a = 3/5 in the 
limit of high Sh p - The growth rate of the m = 2 mode, 

7iin ( m = 2)i nas a c l° se to but somewhat larger than 
1/3. Calculations with m = 4 and m = 6 (not shown 
here) gave intermediate values 1/3 < a m < 3/5 (with 
a m being the scaling exponent for the mode number m 
in the limit of 5h p > S max ). 

In Fig. E| the growth rates of the m = 1 and m = 2 
cigenmodes are compared for Case (lb), where £>i 2 = 
0.21. In contrast to Case (la) [Fig. US D 12 = 0.06], 
now there are two unstable m = 2 eigenmodes. A scal- 
ing exponent ot m =2 ~ 3/5 is obtained for both m = 2 
eigenmodes [while a m= 2 ~ 1/3 in Case (la)]. 

As a result of the complicated dependence of 7n„ 
on the parameter set {m, s\, S2, D12, Sb_ p }, growth rates 
7iin(?7i > 1) m ay rise above "fn n (m = 1) in certain regimes 
of the parameter space. This can be seen in Fig. 1131 

7^(m = 8) > 7 1 ^ ) (m = 1) for 4 x 10 4 < 5 Hp < 2 x 10 8 , 

and 7^(m = 2) > 7 1 ( i 1 n ) (m = 1) for S Hp > 8 x 10 5 (upper 
limit not known). In this regime dispersion relations are 
found to peak at TO pea k > 1- 

Cases with q s = 2 

In Fig. 1151 the Sn p dependence of the linear growth 
rates of q s = 2 DTMs is shown. A case with relatively 
large D12 is plotted in Fig. I15f a). Both m = 2 eigen- 
modes follow the scaling law 7u n oc S^^ rather well. 
When D12 is reduced the growth-rate scalings deviate 
from this power law in a wide range of <Sh p , as can be 



FIG. 15: (Color online) 5*h p dependence of the linear growth 
rates of q s = 2 DTMs. In diagram (a) the growth rates of the 
two m — 2 eigenmodes of Case (Ilia) with D12 = 0.31 are 
shown [cf., Fig. 21a) and (c)]. In (b) the growth rates of the 
unstable m = 2 and m = 8 eigenmodes (both M' 2 ^-type) of 
Case (Illb) are plotted [cf., Fig.[l» and (d)]. 

seen in Fig. I15f b">. While the higher- m mode still ap- 
proaches a m = 3/5, the lower-m mode has a m close to 
1/3. Note the similarity between Fig. 1151 and the q s = 1 
results shown above. 



IV. DISCUSSION AND CONCLUSIONS 

Double tearing modes with high poloidal mode num- 
bers are destabilized when the distance between the res- 
onant surfaces is small. For a given inter-resonance dis- 
tance the mode number of the fastest growing mode 
was observed to shift to lower m when 5h p is increased 
(Fig. llOfl . This is related to the fact that modes with dif- 
ferent mode numbers m approach scaling laws 7i; n oc 
with different exponents a = a™. Moreover, and in 
agreement with earlier works 0,111] , it was found that the 
scaling exponent a m is a function of the inter-resonancc 
distance D i2 (e.g., Figs. Upland ITlJl . 

Linear tearing mode theory predicts two characteristic 
values for the SH P -scaling exponent: a = 1/3 for m = 1 
modes |25| and DTMs on nearby resonant surfaces 0], 
and a = 3/5 for STMs [26|. When the distance D 12 
is increased, DTMs are transformed into STMs and a m 
gradually increases from 1/3 to 3/5 (except for m = 1 
modes) (28|. In the present work, it was observed that 
the transition from a m = 1/3 to 3/5 also occurs when m 
is increased (Figs. IT51 1141 and [T5|l . However, the variation 
of m does not change the character of the DTM mode 
structures, which is in contrast with the above-mentioned 
effect of increasing D12 (Figs. El and 0J. 

For practical reasons, the knowledge of the mode num- 
ber of the dominant mode m pea k in the spectrum 7ii n (w) 
is important, since it determines the structure of the mag- 
netic islands in the early nonlinear regime . The obser- 
vations that TO pca k varies with S*h p and that the growth 
rates "fn n (jn) may approach their respective characteristic 
scaling exponent a m only at very high values of Sh p [cf., 
Figs. 1131 1141 and rrST b)] have the following important con- 
sequence. Results obtained from nonlinear simulations 
run in the collisional regime (e.g., Sh p ~ 10 6 ) for a given 
configuration may not be easily extrapolated to higher- 
Sh p , since magnetic islands with different poloidal mode 
numbers are expected for different Sn p . Note, however, 
that the effective magnetic Reynolds number in the re- 
connection regions may be smaller due to the action of 
micro-turbulence (anomalous resistivity) [39L l40j . 

In summary, the linear instability characteristics of 
DTMs with high poloidal mode numbers m were stud- 
ied numerically. High-m, tearing modes become unstable 
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when two or more resonant surfaces q(r s i) = q{r S 2) = q a 
are formed in a tokamak plasma and when the distance 
between these resonances D12 = \r S 2 — r s i| is still small. 
It was shown that despite the low magnetic shear in the 
vicinity of <7 m i n , modes with high m have high growth 
rates due to the destabilizing effect of small D12. The 
width of the DTM spectrum and mode number of the 
dominant mode were found to be independent of the q- 
value. Broad spectra of unstable DTMs, with dominant 
modes having m > 1, were found in a wide range of mag- 
netic Reynolds numbers, including the regimes in which 
tokamaks operate. Let us note that this result may also 
be applied to configurations with more than two resonant 
surfaces Q and configurations with low magnetic shear 

The findings of this linear study motivate a nonlin- 
ear investigation of DTMs. Nonlinear simulations were 
performed in the past for relatively large inter-resonance 
distances (e.g., [3, IHEl E3)- According to our re- 
sults, during the stage where the inter-resonance distance 
is still small, fast growing high-m DTMs may signifi- 
cantly modify the q profile near (7 m i n and thereby affect 
the long-term evolution Q. This might have important 
implications for the understanding of the sawtooth crash 



D3 D3 LIJ and other applications of DTM dynamics 
mentioned in the introduction, such as the anomalous 
current penetration during the current-ramp up phase, 
off-axis sawteeth or the formation of ITBs. 

An effect that was neglected here but may decouple 
DTMs in practice is differential rotation 0, Isil 
However, its influence decreases with decreasing inter- 
resonance distance. A factor that is expected to be im- 
portant in regions with low magnetic shear near resonant 
surfaces (as in the vicinity of (? m in) is the pressure gra- 
dient ^5|. This and other extensions are left for future 
study. 
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